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A  NEW  METHOD  OP  EVALUATING  DYNAMIC  RESPONSE  OP 
COUNTER -FLOW  AND  PARALLEL-FLOW  HEAT  EXCHANGERS 


By  Henry  M.  Paynter1  and  Yasundo  Takahashi2 
SYNOPSIS 


From  the  exact  solutions  for  the  frequency  response  of  counter-flow 
and  parallel-flow  heat  exchangers,  successive  parameters  are  calculated 
which  give  direct  information  for  the  heat  exchangers  regarding  transient 
responses  as  well  as  frequency  responses.  The  numerical  evaluations  of 
the  parameters  from  the  design  data  of  heat  exchangers  are  generally  very 
simple,  although  the  formulae  themselves  appear  somewhat  involved.  Good 
coincidence  with  measured  transient  responses  is  demonstrated  on  an  example. 

INTRODUCTION 
3  L 

One  of  the  authors  has  published'^*’  analytical  solutions  of  heat 
exchangers.  But  numerical  evaluations  of  these  results  were  not  simple, 
especially  for  tubular  neat  exchangers  (Fig.  1)  because  they  involved  , 
distributed  parameter  systems.  A  new  method  developed  by  the  other  author5 
can  be  applied  to  these  cases  to  obtain  a  numerical  basis  for  dynamic  response 
calculations.  Thus,  for  example,  an  estimation  of  the  transient  response, 
which  otherwise  would  have  required  complicated  calculations0,  can  be  made 
very  easily  from  the  design  parameters  listed  below. 


1  Assistant  Professor  of  Mechanical  Engineering,  Massachusetts  Institute 
of  Technology,  Cambridge,  Massachusetts. 

2  Visiting  Fellow,  Electrical  Engineering,  Massachusetts  Institute  of 
Technology,  Cambridge,  Massachusetts. 

3  "Transfer  Function  Analysis  of  Heat  Exchangers,"  by  Y.  Takahashi  in 
Automatic  and  Manual  Control,  edited  by  A.  Tustin,  Butterworth  Sci. 

Pub.,  1952,  p.  235. 

4  "Regeluechnizche  Eigenschaften  dar  Glerch  -  und  Gegenstromwaraeaus- 
tauschern,"  by  Y.  Takahashi.  Regelungstechnik  Heft  2,  1  Jarg.,  1953,  P»32. 

5  "On  An  Analogy  Between  Stochastic  Processes  and  Certain  Dynamic  Systems," 
by  H.  M.  Paynter.  (Forthcoming  ASME  Paper). 

6  Conduction  of  Heat  in  Solids,  by  H.  S.  Carslaw  and  J.  C.  Jaeger,  Oxford, 
1950,  pp.  325-330. 
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NOMENCLATURE 


A  =  surface  area  of  tube  walls  (fta) 

ai  3  ^  a3  *  ^  (dl) ,  when  both  are  equal,  a  =  ax  =  aa 

V2ca 


Ot,  A 
*x  V.c 


iui 


i-  °^Aa 
1  V3ca 


a  3  **8^8 
Vac2 


(dl) 


B  3  intermediate  parameter  (dl) 


b  =  ♦  ba,  bj. 


J&Aa.  b Vs 

C,  V,  “2  C.  V-,  ^  7T~~ 

n  h  1  C  v 


(dl) 


C  s  tube  or  shell  heat  capacity  per  unit  length  along  the  flow  (Btu/ft.deg.F) 


c  =  specific  heat  of  fluid  (Btu/lb.  deg.F) 

D  3  intermediate  parameter  (dl) 

E  3  ditto  (dl) 

F  3  ditto  (dl) 

f  3  ditto  (dl) 

G  3  transfer  functions 
g  3  intermediate  parameter  (dl) 

H  3  total  length  of  flow  distance  in  the  heat  exchanger  (ft) 

h  3  running  length  along  tube  side  fluid  (ft) 

K  3  intermediate  parameter  (dl) 

k  3  overall  coefficient  of  heat  transfer  (Btu/ft3.min.deg.F) 
L  3  h/v  3  distance-velocity  lag  of  fluids  (min) 

M  3  intermediate  parameter  (dl) 


n  3  numbers  of  lags  (dl) 


r  3  ViAa  (dl) 

8  3  complex  variable  of  Laplace  Transformation  (dl) 


-  3 


T&  3  skew  time  of  step  response  relative  to  Lj,  (dl) 

Td  3  dead  time  of  lag-delay  model  relative  to  Lx  (dl) 
s  time  constant  of  lag  model  relative  to  Lx  (dl) 

Tffl  3  mean  delay  of  step  response  relative  to  Lx  (dl) 

Tr  3  time  constant  of  root  lag  model  relative  to  Lx  (dl) 

Tg  3  dispersion  time  of  step  response  relative  to  Li(dl) 
t  3  running  time  (min) 
v  3  fluid  velocity  (ft/min) 

¥  =  flow  rate  (lb/min) 
x  =  h/H  (dl) 

C7C  s  film  coefficient  of  heat  transfer  (Btu/fta  min.deg.F) 

~  coefficient  of  skew  (dl) 

P  s  (ax  +  aa)/2 

cf  =  In  (Steady  state  change  in  output/steady  state  change  in  input) 
£  ^ax  ”  aa)/2 

Q  -  temperature  of  fluid  (deg.F) 

/*  -  coefficient  of  variance  (dl) 

S  =  a/U  +  a) 

X=  VL i  (dl) 

0  3  pipe  temperature (deg.F) 

Co  s  circular  frequency 

Subscripts  : 

1  3  tube  side 

2  3  shell  side 
h  3  tube 
s  3  shell 


BASIC  ASSUMPTIONS 


1)  System  parameters  are  uniform  and  constant. 

2)  Complete  mixing  in  crossvise  directions  of  each  flow. 

3)  The  heat  conductivities  of  vails  are  either  infinite  in  directions  at  right 
angles  to  the  flov  or  alternatively  assumed  to  be  included  in  film  coefficients, 
and  zero  in  flov  directions. 

4)  There  are  no  internal  sources  or  sinks  of  heat. 

5)  Pure  counter-  and  parallel-flows  (fig.  1)  are  considered. 

FUNDAMENTAL  EQUATIONS 

The  system  parameters  necessary  for  dynamic  response  analysis  under  the 
stated  assumptions  are  the  following  fifteen  (see  also  fig.  1)  t 

Flow  rates  of  fluids  =  Vx  ,  Va 
Specific  heats  of  fluids  -  c*  ,  ca 
Surface  areas  -  Ax  ,  Aa  ,  A 

© 

Film  coefficients  -  °<x  , 

Film  velocities  =  vx,  va 

Flowing  distance  -  H 

Solid  heat  capacities  =  C,  ,  C 

n  © 

These  are  conveniently  grouped  into  the  dimensionless  forms  (defined  above); 

•i*  >  aa*  >  *  bx  ,  ba  ,  bg  ,  r 

Sour  of  them  are  also  conveniently  grouped  into  the  following  dimensionless 
forms'  for  d-c  gain  calculations  and  other  purposes; 

7.  These  cure  also  given  in  the  following  forms: 

_  kA  .  _  kA 

>  a*-«S7 

where  k  (Btu/mn-fta-F)  is  the  overall  coefficient  of  heat  transmission  of 
the  heating  surface.  This  form  can  be  introduced  by  means  of  the  well-known 
law  of  heat  transmission,  which  under  the  assumption  stated  above,  is  written 
as: 

l/kA  =  1/^Ax  +  l/*aAa 

where  l/kA  is  the  equivalent  resistance  to  heat  transmission 


ftiilb 

b 


«  a 


i 


where  b  -  >  bg.  Almost  all  these  parameters  have  been  necessary  for  the 
conventional  steady  state  design  of  heat  exchangers;  for  example,  for  mean 
temperature  difference  calculations. 


The  running  time  t(min)  and  running  distance  along  the  tube-side  fluid 
h(ft)  are  also  expressed  in  the  following  dimensionless  forme; 

^  =  tAi  ,  x  =  h/H 

How  the  simultaneous  equations  to  be  sol  red  ares 

<*„  -  #1) 


-jf  =  bi  (©!  -  J*h)  ♦  ba  (Oa  -  J*h) 

-  *»> *,«,  - «.) 


a) 


9*. 


at  *  ».(®a  -  0.) 


In  these  equations,  ©j.  and  ©g  are  tube-side  a  r>-  shell-side  fluid  temperatures, 
el.  and  0  are  tube  and  shell  temperatures,  the  double  symbol  +  is  -  for  counter 
flow  (Pig.  2)  and  +•  for  parallel  flow  (Fig.  3). 


In  the  following  treatments,  the  Laplace  transform  solutions  of  Equation 
(1)  are  expanded  in  the  following  form1 


G(e) 


yvr_s  + 


t  8 

B 


sa  - 


a 

b 


a3* 


(2) 


where  the  parameters,  cf  ,  T  ,  T  ,  and  T  are  given  in  terms  of  system  constants 
listed  above.  The  symbol  s  is  tne  complex  variable  of  the  Laplace  transformation. 
The  value  and  significance  of  this  representation  has  been  indicated  elsewhere  (5). 
However,  one  may  say  in  summary  that  /measures  the  steady-state  amplitude  ratio 
between  response  and  disturbance,  Tm  measures  the  uean  time  delay  between  response 
and  disturbance,  Ts  defines  the  dispersion  or  attenuation  and  Ta  the  assymmetry 
or  phase  non-linearity.  This  characterisation  is  very  efficient  for  any  physical 
process,  such  as  those  treated  here,  where  the  step  response  is  mono tonic  non- 
decreasing  in  time. 


J 
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COUNTER-FLOW 

The  Laplace  transfom  solution  (transfer  function)  of  Equation  (1)  is: 


f x  +  fa  V  Ux  +  fa)a  -  4gi«a  eoth 


'(fx  ♦  *a)a  -  4Exga 


where 


f  *  * 


f*  - 


ax'  (ba+s)  + 
b  ¥  8 

aa»  (bj,  +  s) 
b  +  s 


-*•  s  r  +  e 


b  +  8 
6 


-  aa’W 
b  +  s 


glf9  -  gx  when  G(s)  is  defined  as 

/Outlet  tempera turew  /Inlet  temperature  v 
'of  tube  side  fluid''  'of  shell  side  fluid' 

Eia  s  Ea  ?or  th®  G(s)  °? 

t 

/Outlet  temperature  %  /  /Inlet  temperature  » 
'of  shell  side  fluid'  'of  tube  side  fluid' 


Case  1 
(Fig. 2) 


Case  2 
(Fig.  2) 


Now,  the  parameters  of  Equation  (2)  are  determined  by  expanding  (2)  and  (3)  in 
series  in  s  and  comparing  the  corresponding  terms.  The  expansion  is  easier  for 
G~X(S)  than  G(S) .  Th»  results  yield  a  solution  for  the  new  parameters  in  the 
symbolic  form; 

*  3  f«  »  aa*  .  V  ,  *1  ,  ba  ,  bg  ,  r) 

Tm  3  <ai’  »  aa’  »  V  »  »  ba  >  bs  »  r) 

Ts  3  f»  (H*  »  aa’  >  a6'  »  ,  ba  ,  bg  ,  r) 

s  ^3  (ax*  >  aa*  »  aB'  >  ^x  t  ^a  »  »  r) 
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Details  of  the  algebraic  reduction  procedure  are  given  in  Appendix  1.  The  most 
simple  relation  is  given  when  solid  capacities  are  neglected,  and  when  the  Vc 
values  and  the  velocities  are  equal  for  both  fluids.  For  this  case,  taking 
J°~  a/(l+a),  where  a  =  kA  as  abscissa,  the  new  parameters  T^  ,  Ts  ,  are 

Wc  * 

plotted  in  terms  of  conventional  relative  statistical  measures,  in  which  we  define 
coefficients  in  the  form 


Coefficient  of  variance: 


Coefficient  of  skew: 


(5) 


A  zero  value  of  /*  means  that  the  time  distribution  has  no  dispersion  about  the 
mean  T^j  a  zero  value  of  <x  signifies  that  the  distribution  is  symmetric  about 
the  mean  Tm.  From  the  plot  we  can  observe  directly  that  when  fi  -  0,  which  occurs 
for  small  sizes,  with  low  overall  efficiencies,  the  time  distribution  for  a  step 
disturbance  is  symmetric  (  oc  =  0)  and  has  the  quickest  response  (minimum  values 
of  Tm  and  f*  ) .  As  /°  increases,  Tm,  Z4  ,  and  oc ,  all  increase  with  /*  and  c*.  , 
becoming  infinitely  large  as  the  length  of  the  exchanger  becomes  infinite. 


PARALLEL  FLOW 


The  Laplace  transform  solution  of  Equation  (1)  is: 

0(a)  *  81>a  «  ~  ■*  g—  £  _ 2  ]  ’*  *ei8a  ] (6) 

■Jit ]  -  f,)»  +  4«,& 

The  symbols  fx>  8i»  8a  aQd  &x»a  are  the  same  as  defined  on  counter-flow, 
see  also  Fig.  3. 

From  this  equation  the  parameters  of  equation  (2)  are  determined  in  the  same 
symbolic  form  as  equation  (4).  Details,  of  algebraic  reduction  procedure  and 
typical  special  cases  are  given  in  Appendix  2. 

NUMERICAL  EXAMPLE 

As  an  example  of  the  application  of  the  formulas  above  to  engineering 
practice,  one  can  consider  the  special  instances  of  counter-flow  and  parallel- 
flow  exchangers  with  the  following  assumed  characteristics: 

ai  -  1.5  ax»  *  6  b  -  27  a  =  4 

s 

aa  =1.5  a8'  =  2  r  =  3  b  =  3 

s 


STEP  RESPONSE  IN  TUBE  INLET  TEMPERATURE 
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Counter  Flow:  Exchanger 

Applying  the  above  values  to  Equations  (9) ,  (12),  and  (13)  we  obtain 

e  *  =  1.65 

T  a  2.28 
a 

T  =  2.16 

O 

T  -  3.08 
& 

The  measure  of  spread  T  and  the  as symmetry  measure  T  can,  as  before,  be 
expressed  In  terms  of  dimensionless  coefficients,  usually* 

Variation  /*  =  T =*  0.948 

Skew  /T  )3  -  2.90 

These  values  can  be  compared,  for  example  to  those  for  a  unit  lag,  whose  transform 
has  the  form: 


(Lag) 


G_(s)  *  _ 1 


1  +  T^s 


where  A  -  1  and  °t  -  2.  The  variation  can  be  matched  by  adding  a  suitable 
time  delay  term,  since  the  transform  for  this  case  becomes 


(Lag  +  Delay) 


Gb(s)  -  _e 


-Tds 


1  + 


for  which 


T.* 


v^T< 


giving  /*■  =  T  Jt*  *  T/  /Td  +  14 
0C  3  Ta3/T83  3  2 

ft  c 

Thus  a  lag-delay  model  with  /*■  =  0.95  and  cC  *  2  would  have  the  fora  shown  in 
Fig.  5. 

However,  the  value  of  oc.  for  the  heat  exchanger  indicates  a  curve  more 
skewed  than  that  expressed  by  a  unit  lag.  Such  a  function  is  found  in  what  toe 
statisticians  would  call  a  "chi-square"  distribution  of  one  degree  of  freedom 
with  the  transform: 


(Root-Lag) 


?*m**m^m ppp^pppwpp 
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having  the  parameters 

*■ 3  V2 

T.  *  V^2 

/*3  J*  -  1.4U 

9* 

It 

n* 

oc  3  2J2  3  2.838 

If  this  root-lag  function  is  delayed  in  addition,  with  the  transform: 

— T,8 

(Root-lag-delay)  G.(s)  ~  e  ^  /  V 1  +•  T„s 


and  parameters 


T  a  T  +  —  T 
a  Ad  2  r 


T 

Ts  a  ?r/j2  /*  =  yi/(l  +  2 

X* 

0C-  2  Vz 


Ta3Tr 


Thus  a  distribution  curve  of  this  form  with  /*-  0.95  and  oC  -  2.828  is  also 
sketched  in  Fig.  5* 


There  are  many  other  possible  distributions,  all  with  the  same  values 
Tm  ,  T-  $  T«  »  but  differing  in  higher  order  terms.  These  will  all,  in 
general,  give  reasonable  approximations  to  the  dynamic  response  of  the  given 
system,  and  many  are  susceptible  to  ready  calculation.  In  the  present 
instance,  those  shown  in  Fig.  5  come  directly  from  functions  which  are  readily 
available  in  tabular  form  and  also  easily  realized  in  computing  networks,  etc. 


Parallel  Flow  Exchanger 

The  general  formulas  (18)  to  (21)  give 
e^  -  0.475 


T  a  2.81 
m 

T  =  1.10 

s 

T  =  1.02 

a 


with  the  coefficients 


3  T  /  T  =  1.10/2.81  =  0.391 

8  IB 

^  3  <Ta  /  V3  3  (1’02/  1*1°)3  3  °»795 
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Now,  it  is  readily  shown  (3)  that  an  n-lag  cascade  (  a  chi-square 
distribution  of  2n  degrees  of  freedom)  has  the  coefficients 

-  l/-/n  ,  OC  =  2/yn 

In  the  special  case  of  n  =  6.5,  with  a  transfer  characteristic  of  the 

fora: 


Gd  (s)  =  1 _ 

(1  +  T^s)6*5 

The  coefficients  become 

/*  -  0.392,  oC=  0.784 
and  the  mean  time  Tm  is  given  by 

T.  =  h 

Accordingly,  a  reasonable  approximation  to  the  parallel  flow  step  response 
characteristic  can  be  found  In  the  chi-square  distribution  of  13  degrees  of 
freedom  with  =  2.81/6.5  =  0.432.  This  is  indicated  in  Fig.  6  and  compared 
with  a  lag-delay  approximation  with  /*■  -  0.391  and  CL~2. 


FREQUENCY  RESPONSE  RESULTS 

In  terms  of  the  same  representation  used  for  transient  response,  the 
frequency  responses  of  heat  exchangers  may  be  found  directly. 


Thus,  if 


G(s)  -  e 


J  -T_s4rr  T_as2  -  T  T_3s3  »-... 


m 


s 


-  e 


S  *  Vs*** 


Even 


-T  s  -  |  T  3s3... 
m  o  a 

"Odd" 


Then,  with  s  -  jo  1 

Amplitude  |g(w)|  =  e  ^  s 


J  *■  2  ^ca<wa  +*  *  • 


Phase  Z  G(p)  =  -  T^o 


+  1  T 

Z  a 


3^3 


so  that  the  set  of  parameters  ^  ,  T  ,T  ,T  ,  ...  describe  the  frequency 
characteristics  simply  and  uniquely.0  Howfver,aspecifi cation  of  only  the  first 
few  parameters  merely  defines  the  low  frequency  behavior,  and  as  the  frequency 
is  increased,  more  of  the  time  constants  will  be  required  to  characterize 
the  behavior. 
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However,  one  can  proceed  as  was  done  in  estimating  the  step  response, 
by  picking  a  suitable  model  using  the  low  frequency  constants  alone,  and  thus 
extrapolate  the  response  to  high  frequencies  under  the  tacit  assumption  that 
the  model  so  chosen  will  behave  at  least  roughly  like  the  prototype  at  higher 
frequencies. 

Then  for  the  counterflow  exchanger,  there  has  been  plotted  in  Pig.  7  the 
predicted  frequency  response  characteristics  for  the  two  previously  determined 
models,  namely 


COUNTER-FLOV  EXCHANGER 


Lag-Delay  Model 

Root-Lag-Delay  Model 

Amplitude 

1//1  +  Vo> 

3  1/ A  +  4. 66  co a 

Ratio 

X/  M  +  T  acoaQ  2 
-  1/(1  +  9.&J1)0’25 

Phase 

tan”1  co 

-  0.12  tv  +  tan'1  2.16a; 

Tdw  +  |  tan  Tr  oj 

31  0.75^4*0.5  tan“13.07o-> 

In  a  directly  similar  fashion,  the  response  characteristics  for  the 
parallel  flow  exchanger  have  also  been  plotted  in  Fig.  8  from  the  formulas 
below. 


PARALLEL-FLOW  EXCHANGER 


Lag-Delay  Model 

Multi-Lag  Model 

Amplitude 

Ratio 

i/Vi  +  V*3 

1 

=  l//l  +  1.21  OJ* 

1/(1  +  T  a  *,3) 13/4 
r 

=  1/(1  +  0.19c  *)3*25 

Phase 

tan  _1 

6.5  tan  -1  (Tro>) 

—  1.71  uj  +■  tan  1  l.lOco 

—  6.5  tan  1  (o.43 c) 

EXPERIMENTAL  CONFIRMATION 

Experiments  under  carefully  controlled  conditions  have  been  made  previously 
by  one  of  the  authors  (1,2)  upon  a  heat  exchanger  model  used  both  in  counter  flow 
and  parallel  flow.  These  yielded  among  other  results  the  response  in  the  shell 
stream  outlet  temperature  to  disturbances  in  the  tube  stream  inlet  temperature. 
These  disturbances  involved  both  stepwise  and  sinusoidal  changes. 
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Results  of  some  of  these  model  tests  are  Indicated  in  Fig.  9  to  12. 

It  is  important  to  stress  that  the  step  responses  and  the  frequency  responses 
represent  data  from  independent  test  procedures  and  therefore  represent,  in  a 
certain  sense  at  least,  independent  physical  data. 

From  blueprint  data  and  direct  measurements  of  the  model  the  basic 
physical  constants  were  obtained.  These  correspond  precisely  to  the  data 
assumed  in  the  numerical  examples  of  the  previous  paragraphs.  However,  the 
surface  conductance  constants  0Ci  and  0C3  were  back-figured,  at  least  in  part, 
from  the  calculated  steady-state  (aero  frequency)  temperature  ratios.  Moreover, 
the  distance  velocity  lag  was  estimated  from  measurements  only  with  tolerable 
accuracy  at  Lj  -  0.6  minutes  with  a  probable  error  of  at  least  0.05  minutes. 

With  these  restrictions  understood,  the  predicted  and  measured  step  and 
frequency  responses  are  depicted  in  Fig.  9  to  12. 
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APPENDIX  1. 


For  counter-flow,  the  parameters  of  equation  (2)  are  given  by; 
<5  3  -1qD0 

\  *  ®lA> 

Ts  =  /( Di/D0)a  -  2Di/D0' 


a 


=  ^ 7  t)3  -t>t>  +  O 


(7) 


etc* 

The  D0  ,  Dj.  ,  Da  ,  D3  are; 

D»  =  ^Ta  lK“  "  B») 


Dx  * 


.  1 


ai»a 


(Mx  +  Bi  *  ) 


(Ma  Ba  ^  |i) 

S^a  (M3  +  B3  *  T-> 


(8) 


etc. 


where  a1>a  -  ax  for  case  1,  a1>s  =  aa  for  case  2  in  Fig.  2.  The  M  and  I 


^  c  c  J 


M3=i  (>)  (i  -i) 


(9) 


'B8  6 


are: 


etc 


-  u  - 


B0  *  £eoth  £ 

Bx  s  £Ej.  (coth  L  -  £  csech*  £.) 

Ba  3  t[la  (coth  £  -  €.  csech*  )  +  £ex*  (  £  coth  £  -  l)}  ^  ^ 

B3  3  £^  (coth  £  -  £  csech*  £  )  ♦  2  E^Ea  £  csech*  £(  £  cotht-1) 
+  Ex^  E.  *  csech*  £(coth  £  -  £  coth*  £  +■  £ 

etc. 

where  6  =  (ax  -  aa)/2  . 

The  Ex  ,  Ea  ,  in  Equation  (6)  are: 


and 


v  _  Kj_  K  -  -la.  _  IL 

El  “  2K.  '  **  2K0  SKo 


K3  KxKa  Kx' 
3  •£*■—  -  ♦ 


(11) 


=*3  ~  2Ka  "  4K0*  ’  K  3 


K„  * 


Kx  3 
Ka  = 


8  £.* 
b 


(12) 


mi  ^(1  +  r  +  !s)t  [<i  ♦  r  ♦  i») 

o/«  i.  m  a  I  fi.t  A  .  % 


_  2iax 


J^lj*  (  -  2  (ax  .  a3) 

=  -  ♦  6  .%«■!-  <h  ♦  m>  -  a 

tr  b 

*  *  "f  (*x  *  aa)  -  2  b^b-  { (*l'  +  aa')  -  (ax  *  aa)] 

b$  s  V.  # 

♦  -§I  (  (»i  *  «a)  -  (“i*  (1  +  r  *  i?  ) 


-2rr(i*^f) 
a  s 


a. 


Sx 
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Given  the  system  parameters,  ve  can  evaluate  M  and  K,  and  from  K  we  can  find 
E,  hence  B.  Applying  these  M  and  B,  the  required  parameters  5  ,  T  ,  T  ,  and 
Ta  are  found  by  Equation  (7) ,  (8) .  These  procedures  and  relations10 get  limpler 
for  special  cases  as  follows; 


Special  Case  It 

If  £  ?  0,  that  is,  ax  =  aa  ,  the  Equation  (10)  may  be  rewritten  as 


Bo  =1 


Bx  -  Ki/12 

Ba  =  Ka/12  -  V/720 

B3  =  Kyl2  -  K^a/360  >  K13/30240 


(13) 


Special  Case  2» 


If  ve  neglect  the  effects  of  solid  capacities,  the  Equations  (11)  and 
(12)  reduce  tot 


Ex  »  <aa  *  fc)  (1  +  r)/(2£>) 


(U) 


Special  Case  3: 


If  £  ==  0  in  Special  Case  2,  the  Equation  (10)  may  be  directly 

given  by* 

B  =  1  +ea/3 


(15) 


Ba  *  (“£  (j  -  (a*  +  £  )3  "  45”  £  *  *  3l5^3(aa+^)a 

B3  =  ^  8^  (*a  +  £  )  [  ’  “45  +  9^  (&a  +  fc  )a  +  3J5  £  * 


2078 


288225 


(aa  +  £  )a  £  aj 
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Special  Case  4: 


If  £  =  0  in  Special  Case  3#  i.e.,  no  solid  capacities  and  ax  3  aa  *  a, 
the  final  results  are  directly  given  by: 


f  ij»  3  o 

e"  (  -  T  T  a  T  ' *  x  _  n  .  x3  /!_  2aa 

—  JS-SL  ♦  )  -  -  (1  *  r)  (90-^ 


(16) 


Special  Case  5: 


Same  as  Special  Case  4,  and  r 
then  we  have: 


1.  Let  us  denote  e 


a 

1  +•  a 


» 


-  1  ♦ 


I 

3 


r 


+ 


B/° 2 

45  (1  -/>) 


(17) 


2/ 3  .  8^  ,  8Z*3  .  32/*3 

27  15(1  -  77  45(1  -/* )  315(1  -f)* 


The  values  of  these  terns  are  shown  in  Fig.  4 
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APPENDIX  2. 


For  parallel-flow  from  Equation  (6)  the  parameters  of  Equation  (2) 
are  determined  as  follows  : 


ax»a 


1 

*  T  b 


sinh  P  e-<* 

f 


T_  =  M,  +  r  -  tx 


V. 


=  -«a  +  Fa 


(18) 


-f-  =  m3  -  f3  ♦  FA  *  -ffcr  -  V  -  3^ 


where  ax,3  =  ax  for  Case  1,  a1>a  =  aa  for  Case  2  in  Fig.  3,  and  fi  ~  (ax  +■  aa)/2. 

The  H  and  F  in  these  results  are: 

M  -  UlU  +  *!■*  .t.  Pal.  -  L 

"1  2  Oh  2b  h 


S 


m  -  &s  a,1  +  a3'  P 

Ma  "  2b  4  ~  2b4  U* 

s 

H  _ +  SaULSal  _  ji 

3  2b  3  2b3  b3 

8 


(19) 


(20) 


Fx  *  Ex  (  P  coth  p  -1) 

*a  s  ®a  (  /*  «>th  /  -1)  -  Exa  [p  coth  f  -  (1  ♦  ~^)j 

F3  =  E3  (  Z3  coth  /*  -1)  -  2ExEg  fp coth  f  -  (1  +  'll)) 

6  a  /  a  <)  2 
+  Ax3  J  (1*  -g-  )  /coth  p  -  (1  +  -tg-  )J 


The  EXJ  Ea,  ...  in  Equation  (20)  are  defined  in  the  same  form  as  Equation  (11) , 
but  K0,  Ex,  ...  in  it  are  given  in  place  of  Equation  (12)  as  follows: 

*.»♦/** 

Kx  -  2(ax  -  aa)  (1  -  r)  -  2(ax-aa)  -gS.  +.  %  ^(a1-aa)(a1»-aa*)-4^aj 
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"a  =  [(l-r-F>  [*^-f  U  -r-^S)  ♦ 


L».',,r  a,')U  i 'd± 


S3  =  2F 


U  -  r  -  ^  /(•»'  -  »a')  -  (aa  -  aa)| 


♦■fl  (1  -  r  -  ^)(al-aa)  -  t  (ax  -  aa)(ai'  -  aa') 
s  tr  tr 

-4r  (ax*  -  aa‘) 


Special  Case  1: 


No  solid  capacities, 


'  -  -/*_  sinh, 

e  -  6  a^a  p 

T  =  l±i  _ 

«  2 


(  /^coth  — 1) 


(l-Z^csech^ 


T  3  =  ~  aa)3  (X-  r)3  fj.  +  .^L2  -  y*coth/*J  coth/3  + 

a  64/3  5  ^  3 

f !  .  lya)2-.  j  (/>  2  C8eChM  _X]  + 

1  ^  '  32/3 4 


f1  *  ( P "V  -1’ 


Special  Case  2  » 


Same  as  Case  1,  and  r  -  1. 

J  -  a  «“/*  8lDh^ 

®  &1»S®  jg 

T_  -  1  ^ 


IB  =  0,  T#  =  0 
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